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Abstract 

In this paper we analyze scalar gravitational perturbations on a Robertson- Walker background 
in the presence of multiple scalar fields that take values on a (geometrically non-trivial) field man- 
ifold during slow-roll inflation. For this purpose modified and generalized slow-roll functions are 
introduced and their properties examined. These functions make it possible to estimate to what 
extent the gravitational potential decouples from the scalar field perturbations. The correlation 
function of the gravitational potential is calculated in an arbitrary state. We argue that using the 
vacuum state seems a reasonable assumption for those perturbations that can be observed in the 
CMBR. Various aspects are illustrated by examples with multiple scalar fields that take values on 
flat and curved manifolds. 
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1 Introduction 



As has been known for a long time, inflation |7|, |T^] offers a mechanism for the production 
of density perturbations, which are supposed to be the seeds for the formation of large 
scale structures in the universe. This mechanism is the magnification of microscopic quan- 
tum fluctuations in the scalar fields present during the inflationary epoch into macroscopic 
matter and metric perturbations. Also, since a part of the primordial spectrum of den- 
sity perturbations is observed in the cosmic microwave background radiation (CMBR), this 
mechanism offers one of the most important ways of checking and constraining possible 
models of inflation, see e.g. ||. 

The theory of the production of density perturbations in the case of a single real scalar 



field has been studied for a long time 19, H|. However, to realize inflation that leads to the 
observed density perturbations in a model without very unnatural values of the parameters, 
it is now thought that one needs more than one field. This is a strong motivation for hybrid 
inflation models fl3| (more models can be found in [lj]). Also, many high-energy theories 
contain a lot of scalar fields. The Higgs sector of the standard model consists of one physical 
particle, but in grand unification or super symmetric models one expects many more scalars. 
Ultimately one would hope to be able to identify those fields that could act as inflatons. 
For all these reasons it is important to develop a theory for perturbations from multiple 
field inflation as well. Work in this direction has been done by several people. Using 
gauge invariant variables the authors of |^3|, H treated two field inflation. The fluid flow 
approach was extended to multiple fields in |14|| , while a more geometrical approach was 
used in |26|, both methods assumed several slow-roll-like conditions on the potential. 
Using slow-roll approximations for both the background and the perturbation equations the 
authors of [p0| , |24[| were able to find expressions for the metric perturbations in multiple 
field inflation. 

In this paper we generalize the slow-roll parameters for a single background field to 
the multiple field case in a systematic way. We can then give a clear quantification of the 
relative importance of terms in the equations obtained by extending the single field density 
perturbation calculations by Mukhanov, Feldman and Brandenberger to multiple fields. 
In our definition of the slow-roll functions we do not implicitly assume slow roll to be valid. 
As a consequence, these slow-roll functions are expressed in terms of derivatives of the 
field velocity and the Hubble parameter, but not the potential as the conventional slow-roll 
parameters are. A big advantage of this is that these slow-roll functions can be identified 
in all kinds of equations that are still exact, which is the reason why we can estimate the 
relative importance of various terms and make a well-motivated decision about neglecting 
some of them. In the case of multiple scalar field inflation it is very convenient to think in 
terms of vectors: if the fields are local coordinates on a curved manifold, their derivatives 
and their fluctuations can be interpreted as covariant vectors in the tangent bundle of the 
manifold. Since the fluctuations are assumed to be small, a linearization procedure can 
be used to obtain equations for the perturbations. In these equations a prominent role is 
played by the scalar field velocity. Since in single scalar field inflation the only direction is 
parallel to the field velocity, it is to be expected that the parallel field perturbations can 
be absorbed in the gravitational perturbations, as happens in the single field theory. Using 
the modified definitions of the slow-roll functions we can investigate to what extent this is 
correct. 

In the CMBR spectrum we can observe correlations in the temperature distribution. 
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They are assumed to be due to gravitational perturbations that are of a quantum origin 
at the beginning of inflation. To calculate these correlations we therefore need to address 
the question of what is the relevant state at the initial stages of inflation. We argue that 
the conventionally used vacuum state seems a good assumption. We investigate the effect 
on the correlator in the (probably unrealistic) case that the true state at the beginning 
of inflation is thermal with a temperature of the order of the Planck scale. For most of 
this work we restrict ourselves to the calculation of the correlation function near the end 
of inflation. However, to show that our results are consistent with other results in the 
literature, the correlation function is also evaluated at the time of recombination, ignoring 
possible problems during (p)reheating. 

During inflation there is a relatively sharp transition in the behaviour of a fluctuation 
when the corresponding wavelength 'passes through the horizon' (this is discussed in detail 
in this paper). This moment of horizon crossing can be used to identify a certain scale k. 
The smallest scale (or largest wavelength) that can be observed in the CMBR is the one 
that is reentering the horizon at this very moment, indicated by ko. Assuming that the 
fluctuations do not change once they have passed outside the horizon,^ we can observe the 
undisturbed inflationary perturbation spectrum up to the scale that reentered the horizon 
at the time of recombination when the CMBR was formed, which has ln(/c rec /&o) ~ 3.5. 
On the other hand, larger scales already reentered the horizon before this time, and so 
on these scales the spectrum has been influenced by physical processes taking place long 
after inflation. The largest scale expected to be measured with the Planck satellite has 
ln(/c/feo) ~ 7 (corresponding with multipole / = 2000). An important quantity in our 
discussions is the number of e-folds Nk that occur after a certain scale k crosses the horizon 
during inflation, until the end of inflation. As is derived in e.g. equation (5.16) of fpTOf] , 
depends logarithmically on k, taking a value of about 60 for k^. This number has only 
a logarithmic dependence on model-dependent quantities like the reheating temperature. 
Hence the observationally important scales cross the horizon about 50 to 60 e-folds before 
the end of inflation. 

Apart from this introduction and the conclusions at the end, the paper is structured 
as follows. In section [2] the theory of scalar perturbations with a gravitational potential 
coupled to a single scalar field is reviewed, following the methods discussed in |j9j| . This 
allows us to introduce various relevant concepts. In addition, this section also makes it 
possible to compare the single and multiple scalar field situations in the following sections. 

The generalization to multiple scalar field inflation is developed in section |3[ As these 
scalar fields parameterize a possibly curved manifold, it is necessary and convenient to in- 



troduce some geometrical tools in subsection 3.1. Two examples of these are the inner 



product and the covariant derivatives associated with the metric of the manifold. Another 
concept introduced here is the projection on directions parallel and perpendicular to the 
background field velocity, which plays an important role in the decoupling of the gravita- 
tional perturbations from the independent scalar field fluctuations. In the next subsection 
we explain why the single field method does not lead to an immediate decoupling of the 
gravitational perturbations in the multiple field case. To be able to compare the relative 
size of terms in the equations we define modified versions of the slow-roll functions that 



1 In ]29| it was proved that this is always true for the so-called adiabatic perturbations, but it need not 
be true for the so-called isocurvature perturbations, see Q . A full treatment of this complication is beyond 
the scope of this paper. 
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can be used in multiple field inflation in subsection 3J and derive various useful properties 
for them. In subsection [O] we derive the equation of motion for the gravitational potential 
in terms of these slow-roll functions and show that the term coupling this equation to the 
perpendicular field perturbations is small, so that also in the multiple field case there is an 
effective decoupling. The equation of motion for the perpendicular scalar field perturba- 
tions is derived in the appendix using similar methods. A discussion of the solution for the 
decoupled gravitational potential concludes this subsection. 

Section || is devoted to the computation of the quantum correlation function of the 
gravitational potential. We argue why taking a vacuum state at the beginning of inflation 
to evaluate this correlator seems a good approximation. 

In section |B| we analyze the behaviour of the background scalar fields during slow-roll 
inflation in various multiple field cases. In subsection 5T two examples of a quadratic 
potential on a flat manifold are considered: with equal and with different masses. Next, we 
turn to the generalized situation of a curved manifold with arbitrary potential in subsection 
|5.2| . To illustrate some aspects we take the curved manifold to be the sphere with embedding 
coordinates. In each of these examples we compute the model dependent factor that appears 
in the expression for the correlation function of the gravitational potential. 
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2 Scalar perturbations 



2.1 Linearized gravitational perturbations 

We can divide general perturbations in the universe in the following two types: metric per- 
turbations and matter perturbations. Of course they are related by the Einstein equations. 
On the other hand we can divide both matter and metric perturbations in three different 
classes: scalar, vector, and tensor perturbations [§, |28|, depending on how they trans- 
form under spatial transformations of the background metric. In this paper we consider 
only scalar perturbations since they are the main cause of the fluctuations in the CMBR. 
We assume throughout this paper that all perturbations are small, as on the one hand they 
presumably originate from quantum perturbations, while on the other hand the fluctua- 
tions in the CMBR that we observe are tiny. In particular this means that we linearize all 
equations with respect to the perturbations. This section is essentially a review of |H| . 

The Robertson- Walker metric for a spatially flat background combined with scalar met- 
ric perturbations may be written as 

ffiiv, x) = 9,M + Sg„(v, x) = a 2 ("q 1 + a 2 ^ ^ J . (1) 

We take a flat background with scale factor a(rj), since we are going to apply our formulae 
to the universe during and after inflation, when it has already been inflated to complete 
flatness. <£, B, E are four scalar functions of spatial coordinates x and conformal time rj 
which together describe the metric perturbations. One often refers to as the Newtonian 
potential, as the 00-component of the metric in a weak field approximation can be identified 
with the potential of Newtonian gravity. The conformal time r/ is related to the comoving 
time t by dt = adrj. The advantage of using conformal time is that then the scale factor a 
is an overall factor of the full metric, not just of the spatial part. We take the scale factor 
to have dimension of length and therefore r\ and x are dimensionless. Differentiation with 
respect to conformal and comoving time are denoted by ' = and ' = dt, respectively. The 
conformal Hubble parameter Ti is defined as Ti = a' /a. It is related to the comoving Hubble 
parameter H = a/ a by Ti = aH. The notation j denotes a derivative with respect to the 
spatial coordinate x l . 

There is a problem with the interpretation of the metric perturbations because it is 
difficult to separate the physical metric perturbations from the ones that can be gauged 
away by a coordinate transformation. Of course the final, physical results do not depend 
on the choice of coordinates. However, it is often convenient if intermediate results are 
also independent of the gauge chosen. To this end so-called gauge-invariant quantities are 
introduced j|, which are defined to be gauge invariant with respect to infinitesimal 
coordinate redefinitions. In this approach one defines gauge- invariant metric and matter 
quantities and uses those in the full system of metric and matter perturbations^] They are: 

$0/0 = $ + ! \(B-E')a]', *k?0 = t!-H(B-E'), 5q i9i) = 5q + q'(B - E'). (2) 
a 

Here we have separated an arbitrary scalar (matter) quantity qmiiVi x ) into a homogeneous 
background part q(rj) and a perturbation 5q(rj,x) = qi u i\(rj, x) — q(rj), just as we did for the 

2 The main alternative is the fluid flow approach jl^, where one makes certain assumptions about the 
matter content of the universe, which enables one to eliminate the metric perturbations and derive equations 
for the matter perturbations in closed form. 
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metric. As can be seen from (g), working with these gauge-invariant quantities is equivalent 
to choosing the longitudinal gauge in which B = E = 0. As we only employ the longitudinal 
gauge in this paper, we omit the (gi) labels without risk of confusion. 

For most of the treatment in this paper it is irrelevant whether scalar perturbations like 
6q(j],x) and <5(7/,x) are classical or quantum objects. This is because we will linearize in 
those quantities so that the quantum nature (such as variables that do not commute) does 
not play a role. Hence we may derive and manipulate the equations as if all quantities were 
classical. Only when we are computing the correlator of the Newtonian potential ($<&) in 
section || do we have to take the quantum nature of the perturbations into account. 

We now treat the dynamics of the universe with scalar perturbations. The background 
Einstein equations read 

H 2 = -I K VT£, (2H' + H 2 ) 5) = -K 2 a 2 T), (3) 

and T 1 - = 0. Expanding the Einstein equations to first order in the perturbations gives 

-3H{H$ + + Air = --K 2 a 2 ST° , (4) 

+ = -^ 2 « 2 <5T°, (5) 



6 j- = l -K 2 a 2 5T) } (6) 



with k 2 = 8itG = 87I-/MJ, and A = Y^i d i d i- Later 

on we often switch to complex Fourier 

modes fk(r]), defined by 

f(v, x) = J d 3 k (Afaje-** + /^y kx ) , (7) 

where / is any real quantity that depends on both time and space coordinates, e.g. $(?7,x). 
After this switch equations for f(rj,x) become equations for /k(??) and the spatial Laplacian 
—A is replaced by k 2 = |k| 2 . 

The complicated system of Einstein equations is simplified considerably when the matter 
is described by a scalar field theory. For a scalar field theory with an arbitrary number of 
fields one can easily verify that <5T*- oc 5j to first order in the perturbations, while 5T® can be 
written as 5T® = 5F^ to the same order. Here 5F is a scalar function of the fields and their 
perturbations which depends on the scalar field theory. A very important simplification 



then follows from an argument by Mukhanov et al. IE], who show that one can take = $ 



ifT^. 



By taking a normalized trace of the (r/)-components and subtracting the (OO)-component 
of the background Einstein equations (||), one obtains the following equation: 

n 2 -H' = \K 2 a 2 ^-T\-T^. (8) 

A similar procedure with the perturbed components of the Einstein equations (f|) and @, 
setting ^ = gives the equation of motion for the Newtonian potential $ 

+ Mi& + 2(H' + 2H 2 )$ - A$ = X -K 2 a 2 QtfTj + 5T° Q ^j . (9) 
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With a suitable choice for the Newtonian potential $ to eliminate the integration constant 
the perturbed (Oz)-component (||) of the Einstein equations can be integrated using 5T® = 

+ H$ = -^K 2 a 2 5F. (10) 

In the next subsection we apply these equations to the case of a scalar field theory with a 
single field, and in section to the multiple field case. 



2.2 Scalar perturbations due to a single scalar field 

Now we consider the case of a single real scalar field, both to review the method of Mukhanov 
et al. [19] and to be able to compare it with the multiple scalar field case when we treat the 
latter in section @. From the Lagrangean 



C = ^[--d^d^-V{<t>)) (11) 



we obtain the equation of motion for the scalar field, 

D^cf> - V+ = 0, (12) 
where is the covariant space-time derivative, and the energy-momentum tensor 

T„ = " K (^d X <t>dx<l> + V{<j>)\ . (13) 

In these three equations, <j> denotes the total field </>f u ii. Now we separate the background 
from the perturbations, as defined below (Q), so that in the remaining equations 4> denotes 
the background part of the field. The integrated (Oz)-Einstein equation (|lG|), 

+ = i/cV<ty, (14) 

and the background equation of motion for the scalar field, 

cf>" + 2H(f>' + a 2 V^ = 0, (15) 

can be used to eliminate the fluctuation 6(f) and the first derivative of the potential V<j> from 
the equation of motion (||) for the Newtonian potential 

+ 6H& + 2(H' + 2H 2 )<$> - A$ = -K 2 a 2 V^4>. (16) 

Then this last equation takes the form of a homogeneous differential equation: 

+ 2 (n - ^-J + 2 \H' - H ^ J $ - A$ = 0. (17) 



Using equation (B), 



n 2 -n' = ^V 2 , (is) 
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equation ( |l7[) can be written in terms of u = -^tti^ as 

nil nj 

u" - -r u - Au = 0, with = —. (19) 

6 aqy 

The factor k -3 in the definition of u has been introduced to give it mass dimension one. 
We can use ([18]) and the relation 7i = aH to obtain expressions for u = &/(k 2 \/ —2H) 
and 6 = nH/(a\f— ~2H) that do not contain the scalar field 0, but only quantities that are 
well-defined also after inflation. Equation (|l^) for the variable u is very important because 
it can be used throughout the evolution of the universe: during scalar field, radiation, 
and matter domination |T^| (only during matter domination an extra factor containing the 
sound velocity has to be added). 

By varying (|j~2| ) to first order in the perturbations we find the equation of motion for 
the scalar field fluctuations: 

(aj + 2Hd v - A + c?V M ) 5<P = 40'$' - 2a 2 V^. (20) 

Notice that we did not need to use this equation in our derivation of a homogeneous equation 
for $ (|l7|) or u (|l9|). It was not needed, because the equation of motion of the scalar field 
can be derived from the constraint that the energy-momentum tensor is diver genceless, 
D p T^ v = 0, and is therefore not an independent equation. This is closely related to the 
fact that we could solve for Scj) by dividing the integrated (Oz)-Einstein equation ( |l4|) by the 
velocity cp'. In the case of more fields this constraint can no longer reproduce all equations 
of motion, so we expect to need the equations of motion for the field perturbations in that 
case. 
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3 Inflation with multiple scalar fields 
3.1 Geometrical concepts 

We now turn to the multiple scalar field case, where the scalars 4> = (<P a ) can be interpreted 
as the coordinates of a real manifold Ai on which a metric G is defined. To make optimal 
use of the geometrical structure of this manifold in our discussion of the dynamics of scalar 
fields and their perturbations, we need to introduce some geometrical concepts. From the 
components of metric G ab the metric-connection T bc is obtained using the metric postulate. 
The definition of the manifold Ai is coordinate independent, therefore the description of 
this manifold is invariant under non-singular local coordinate transformations 

cj) a — ► = X a (4>). (21) 

In our treatment of scalar perturbations due to multiple scalar fields we heavily rely on the 
concept of tangent vectors. A vector A = (A a ) is called a vector in the tangent space T p Ai 
at a point p G Ai if it transforms as 

A* A a = X a b ( ( f } )A b , X a b (4>) = X° b {<f>), (22) 

where the comma denotes differentiation with respect to local coordinates. A simple ex- 
ample of a (tangent) vector is the differential d<j>. The cotangent space is the dual of the 
tangent space. Its elements are linear operators on the tangent space 



*C: T P M 
A 




As C a A a is a scalar object, the cotangent vector *C transforms as 

C a — C a = C b (X^) b a . (24) 

The metric G a b can be used to construct a cotangent vector (AJ) a = A b Gb a from the tangent 
vector A. Using index-free notation this reads At = A T G. The notion of (co)tangent 
vectors defined at a point p 6 Ai can be extended over the whole manifold Ai by interpreting 
them as sections of the (co)tangent bundle. 

Using the metric G we introduce an inner product of two vectors A and B on the 
tangent space of the manifold and the corresponding norm 

A B = A f B = A T GB = A a G ab B b , \A\ = VA • A. (25) 

The Hermitean conjugate Lt of a linear operator L : T p Ai — ► T p Ai with respect to this 
inner product is defined by 

B • (VA) = (LB) • A, (26) 

so that Lt = G _1 L T G. A Hermitean operator H satisfies Ht = H. An important exam- 
ple of Hermitean operators are the projection operators. Apart from being Hermitean, a 
projection operator P is idempotent: P 2 = P. 

To complete our discussion on the geometry of Ai we introduce different types of deriva- 
tives. In the first place we have the covariant derivative on the manifold, denoted by - a or 
V , which acts in the usual way, i.e. 

A a . b = V b A a = Ay + T a bc A c (27) 
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on a vector A a , and V- a = V a y = V A on a scalar V. It is convenient to also introduce 
index-free notation for (covariant) derivatives. On a scalar function V (e.g. the potential), 
the derivative d and the covariant derivative V are equal 

(W)« = (dV) a = V a . (28) 

Since we represent dcf) as a standing vector, V and d are naturally lying vectors and therefore 
V T and d T are standing vectors. The second covariant derivative of a scalar function V is 
a matrix with two lower indices: 

(V T W) ab = V a V b V. (29) 

Of course, the same holds good for ordinary derivatives d. 

The curvature tensor of the manifold can be introduced using tangent vectors B, C, D: 

[R(B, C)Df = R a bcd B b C c D d = (r^ c - T a bc , d + Tt d T a ce - Tt c T a de ) B b C c D d . (30) 

One should realize that for later notational convenience we do not use the standard definition 



as made for example in [21|: our R(B,C)D is conventionally denoted by R(C,D,B). 

Next we discuss how spacetime derivatives act on spacetime dependent tangent vectors 
and their derivatives. Purely spacetime covariant derivatives are denoted by and are 
defined in the usual way. The covariant derivative T>^ on a vector A of the tangent bundle 
is defined in components as 

V„A a = d^A a + T a bc d^ b A c , (31) 

while T>n acting on a scalar is simply equal to d^. 

After the introduction of this standard geometrical machinery, we now develop some 
concepts to describe a time dependent scalar field background. They are used when we 
consider multiple field inflation. Consider a curve 4>{t) on a manifold A4, parameterized by 
a real variable t. In later sections, when we describe the time evolution of coupled systems 
consisting of multiple scalar fields and gravity, this variable t is interpreted as comoving 
time. Along this curve the nth derivative vector can be defined by 

0(!) = = ^ and 0< n ) = V^ 1] 4> for n > 2. (32) 

In applications in later sections cj>^ and <fi^ represent the velocity and acceleration of the 
background scalar fields, respectively. In general the vectors cj>^ , cj)^ , . . . do not point in 
the same direction. From these vectors a set of orthonormal unit vectors is obtained by 
using the Gram-Schmidt orthogonalization process. The first unit vector ei is given by 
the direction of cf>^\ The second unit vector e2 is determined by that part of 0^ 2 ^ that is 
perpendicular to ej, and so on. To obtain the direction of cj)^ perpendicular to ei, we use 
the projection operators Pi and Pj^ that project on subspaces parallel and perpendicular 
to (f>^\ respectively, and require that e2 is proportional to P^cj)( 2 \ 

These definitions can be extended to e n , P n , etc., for any n. The unit vector e n points in 
the direction of cf)^ perpendicular to the first n — 1 unit vectors ei, . . . , e n _i. The operator 
P„ projects on e n and P^; projects on the subspace which is perpendicular to ei, . . . ,e n . 
This subspace is also perpendicular to the derivative vectors <ffi~\ ■ ■ ■ , (f)^ ■ To obtain all 
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these objects at once, we let Pg = 1 be the identity and define the mutually orthonormal 
unit vectors e n (for n = 1, 2, . . . ) from <p( n ' and the projection operator P^-i by 



p-L (k( n ) 

^ = |Pn-i0 (n) l, *n= n 'Z , Pn = e n eL P^ = 1 - £ P ? . (33) 
By construction the vector (j>( n ' can be expanded in these unit vectors as 

n 

= ( Pl + . . . + P n ) = £ ^)e p , = e p • 0W. (34) 

P =i 

In particular, we have that (j>n = e n • cf)^ = |P^_i0^|. As the projection operators Pi 
and P^ will occur frequently in later sections, we introduce the short-hand notation: 

pll=p 1 = M-, p-L = pj- = t -pll. (35) 

|0|2' 1 V > 

In terms of these two operators we can write a general vector and matrix as follows: 

A = A II +A ± , M = M II +M |i +M i| +M li , (36) 

with A" = P"A and M" H = P"MP", etc. Notice that because of the hermiticity of the 
projection operator A 1 - • B = A • B -1 - = A 1 - ■ B -1 -. 



3.2 Multiple scalar fields and gravitational perturbations 

We now consider scalar fields <fi that are the local coordinates of a manifold M, using the 
geometrical concepts introduced in the previous section. The space-time derivative of the 
background field d^cf) transforms as a vector, even though the fields <j> in general do not, as 
they are coordinates on a manifold. Also the field perturbation Scf) and its gauge-invariant 
form defined by (pl) transform as vectors. Only the covariance with respect to the coordinate 
transformations ( pi|) of the target space M is manifest in our treatment, because we use a 
flat Robertson- Walker background and work to first order in the perturbations. 

The Lagrangean for the scalar field theory with potential V on the manifold A4 can be 
written as 

£ M = V=9 (-\&4> ■ W ~ y) • (37) 
The equations of motion for the scalars are given by 

ST - P*,) d x cj> - G 1 V T V = 0, (38) 

and the energy-momentum tensor is 

Tt = d»<f> ■ d v <t> - 8* Q<9 A • d x cj> + Vj . (39) 

In these three equations, <fi denotes the total field cf> {ull . Now we separate the background 
from the perturbations, as defined below (^), so that from now on cf> will always denote the 



10 



background part of the field. In this case of multiple scalar fields the equation of motion 
(|9|) of the Newtonian potential $ reads to first order in the perturbations 

$" + 6H& + 2(W + 2H 2 )§ - A$ = -« 2 a 2 (W 5(f)). (40) 

The integrated (Oz)-component ( |i"0| ) of the Einstein equations takes the form 

$' + H$ = -K 2 <p' -Scf). (41) 

In this case it is not possible to construct the analogue of the single field homogeneous 
equation of motion (17) for the Newtonian potential <£, because (41) no longer contains 
a simple multiplication of two scalars, but an inner product of two vectors. Therefore, 
one cannot extract an explicit expression for S(j), as was the case in equation ([l4|) for 
the single field situation. To overcome this difficulty, we divide the field perturbation Scf) 
in a part that is parallel to the velocity field cf>' and a part that is perpendicular, using 
the projection operators (|35|) defined in the previous section.^ Here we use that for the 
projection operators P" and P -1- there is no difference between using comoving time or 
conformal time: they only depend on the direction of (f), which is the same as the direction 
of 4>'. Once we have separated the fields in this way, the parallel part d<fi^ can be eliminated 
in a way analogous to the single field case. 

Using the integrated (Oi)-component of the Einstein equations (|4l|) together with the 
background equation of motion for the scalar fields, 

V^' + 2H(f)' + a 2 G~ 1 V T V = 0, (42) 



the right-hand side of equation (40) for $ can be rewritten as 
-K 2 a 2 {W8cj>) = -K 2 a 2 (G _1 V T F) • 8<p = k 2 (V v cf)' + 2Htf>') ■ (ei(ei • 6<f>) + Scf)^ 
= 2(*' + W$) L^-AV^') ■ ei + ZHj + K 2 (V v cf>') ■ 5^, (43) 

where we used the definition of the projection operators (|35|). Inserting this expression in 
(|40|) and realizing that \ (f>'\ = [V^cj)') ■ <fi', we get 



+ 2 ( H - Jj^D + 2 (%' - ^Trir) ^ - A$ = k 2 (V v 4>') ■ 8^. (44) 



Notice that, apart from the right-hand side, equation (|44| ) looks identical to equation 
for the single field case. However, it is exactly this right-hand side which makes it necessary 
to solve a coupled system of differential equations, as opposed to the decoupled system 
in the single field case. In this case we need the equation of motion for the scalar field 
fluctuations, 

(V 2 + 2HV ri - A - R(0', cf>') + a 2 M 2 (cf>)) Scj) = 4$>' - 2a 2 $ G" 1 V T V, (45) 
where we have introduced the mass-matrix 

M 2 = G- X V T W, (M 2 ) a b (</>) = G^i&VM. (46) 



3 A similar decomposition in the case of two field inflation was simultaneously developed in 
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Notice that R(0', 0') = a 2 R(0, 4>), so that it is possible to absorb the curvature term 
into an effective mass matrix. Equation ([Iq ) is the multiple field generalization of (|20|), 
However, we do not need the total perturbations, but only the perpendicular part, as can 
be seen from (f44|). This system of equations for <& and 8cf>^~ we will analyze in section 3.4 , 
but before we do that we consider the background equations during slow-roll inflation. 



3.3 Multiple field slow-roll functions 

Slow-roll inflation is driven by a flat scalar field potential that acts as an effective cosmo- 
logical constant because of the small slope. In the case of a single scalar field, the notion 



of slow roll is well-established (see e.g. [10, 14, 11]). In this paper we generalize this con- 
cept to multiple scalar fields in a geometrical way. Afterwards we discuss how our slow-roll 
functions are related to the well-known single field slow-roll parameters. 

To define the slow-roll functions we use comoving time t, since then the background 
equation does not contain the rapidly changing scale factor a. The background equation of 
motion (|42j), the Friedmann equation pf), and equation (H) read in comoving time 



V t 4> + 3H(f> + G 1 V T V 



0, 



H 2 



r 



h<t>\ 2 +v 



H 



\A<t>\ 2 - (47) 



The system is said to be in the slow-roll regime if \T>t4>\ <C |3//<^| and ^\4>\ 2 <C V. A more 
precise definition is given below (^lj). Using ( |47| ) the last condition can also be written as 
(-H) « \H\ 



The functions 



H 

w 



(2) 



ff|0f 



and = 



0g) 
H*\j>\ 



(48) 



can be defined whether or not slow roll is valid. Since both i) and £ are vectors, they can be 



decomposed in components using the unit vectors defined in section 3.1. The components 
of these vectors in the directions ei,e2 and e3 are given by 



e 1 



ei • rj 



Vt4>-ct> 
H\<f>\ 2 ' 

V 2 t<t> ■ 
H 2 \c(>\ 2 ' 



e 2 • rj 



6 = e 2 • £ 



H\cf>\ 

e 2 • V 2 4> 
H*\j>\ ' 



& = e 3 • £ 



(49) 



Since £ in general has two directions perpendicular to ei, we cannot use the ambiguous 
notation £ . However, since £ is a vector, £ is defined. The components of £ are not 
needed for the background equations, but will appear in the equations for the perturbations, 
see section |3.4| . The functions e, fj^, and £3 are non-negative, while fj^, fH, and £2 can also 



be negative (apart from the sign, fj" is equal to 



H\4>\ 



In terms of these functions the Friedmann equation ( f47j) reads 



H 



K 

75 



3 



-1/2 



(50) 
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For a positive potential V the function e < 3, as can be seen from its definition. Inserting 
this expression and the functions fl4S|) into the background equation gives 



+ 



2 1 

y/3 K 



G~ l V T VV 




f] + 



ieei 



(51) 



From this equation and the one above, both of which are still exact, we can define precisely 
what is meant by slow roll. Slow roll is valid if e, y/ifj^ and y/ifj 1 - are (much) smaller than 
unity. For this reason e, fj" and fj 1 - are called slow-roll functions. However, 77" and fj 1 - could 
even be somewhat larger than one during slow roll, if e is sufficiently small. On the other 
hand, we will often use the somewhat stronger condition that e, fj" and fj ± have to be small 
individually. The components of £ are called second order slow-roll functions, and they are 
assumed to be of an order comparable to e 2 , er/H, etc. If slow roll is valid, we can expand in 
powers of these slow-roll functions. For example, we can expand the previous equation to 
lowest non-zero order in slow roll, which gives 



+ 



2 1 

y/3 « 



eei + 7/ + 0(fj"e,fj e, e 



~2\ 



(52) 



Next we derive some useful expressions for the slow-roll functions in terms of conformal 

H(l + T?"), 



time. Using (49) and (|35| ) we find 

I0T 



I0'l 

n' = n 2 (i-e), 



\d>'\ 

I0T 



H 2 (3r) + er, 



(53) 



I0'l 



I0'l 



H [2 — e + 377" + (fj ) + ^ 



Here we differentiated the expression = ■ 0' with respect to r] to find an 

expression for |0'|". Differentiating the slow-roll functions with respect to conformal time 
r\ and using some of the above results we find 

e' = 2We(e + r?ll), (rjW)' = H$ + (r^) 2 + ef?" - ^«) a ], 25,^ = W[| + (e - (54) 

The slow-roll functions ([S]) are all defined as functions of covariant derivatives of the 
velocity 0, the velocity cf> itself, and the Hubble parameter H. If the zeroth order slow-roll 
approximation works well, that is if the right-hand side of ([H]) can be neglected, as well as 
the e in fl50D, then we can use these two equations to eliminate cj> an d H in favour of the 
potential V. This is the way the conventional single field slow-roll parameters are defined. 
However, this conventional definition has the disadvantage that the slow-roll conditions 
become consistency checks. Hence, while we can expand the exact equations in powers of 
the slow-roll functions, that is impossible by construction with the conventional slow-roll 
parameters.^ Therefore, we only show what our slow-roll functions look like in terms of the 

4 In the context of single field inflation this fact was noted before and discussed in detail in |ff|] . 
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potential in this approximation, but we do not adopt it as the definition: 



1 |W| 



1 VVM 2 G _1 V T V 



tr 



fM 2 



K 



v\vv\ 



V 



1 \P ± M Z G- L V 1 V 

~t$ VWv\ 



tr (M 2 ) !I± (M 2 )" 



V 



(55) 



The effective mass matrix M 2 is defined in (^6|). Here P" projects along the direction 
determined by W, which is to lowest order identical to the direction of 4> ■ 

In order to avoid confusion and for later comparison, we finish this section by very 
explicitly comparing the slow-roll functions we defined in Jig ) with the ones conventionally 
used in the single field case, e and rj: 

1 V d> 1 V ,<t><P ~|l , ~ tKR\ 

e= 2^Vl = e > ?? = ^^ = - ??+e ' (56) 

where the last equalities in both equations are only valid to lowest order in the slow-roll 
approximation. Of course, fj 1 - does not exist in the single field case, as there are no other 
directions than the parallel one. 



3.4 Decoupling of the perturbation equations 

In this section we analyze the perturbation equations. Starting point are equations (|44|) for 
$ and (|4^) for 5<fi. First we rewrite these equations in such a way that we can draw some 
important conclusions using the slow-roll functions defined in the previous section. It turns 
out that the equation for the (redefined) gravitational potential decouples from the field 
perturbations up to first order in slow roll. Next we concentrate on solving this equation. 

We can write the system of the perturbation equations as a homogeneous matrix equa- 
tion for the vector 5<j> ). However, to remove the first order derivative term from the 
equation of motion for $ we define new variables: 

u = ^VTn^ 1 ^ = -^——r^ 1 ^^ Sv = -TTTn S( t >± = ~7= "/- S^- (57) 

The additional factor of k _1 in the definition of u has been included to make the mass 
dimension of both u and <5v equal to one. Another important point to note is that we chose 
our redefinitions in such a way that no relative slow-roll factors have been introduced in the 
relation between u and Sv as compared to the relation between $ and <50 _L . 

We derive the u component of the matrix equation here, as it is the component that 
we use most. The equation for Sv is more complicated, which is why we only give the 
result here, and refer the reader to the appendix for the derivation. We start with rewriting 
equation ( [44[) in terms of u and obtain 

Using the conformal time version of the third equation in (|47|), 

H 2 -H' = ^K 2 \cf)'\ 2 , (59) 
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and the definitions of Sv and the slow-roll functions (|5^) we can rewrite this as 

u" -Au + H 2 (-21 - ?]H - 2(r/H) 2 + (fi L f + £H) u = V^H 2 ^^ 1 - • <5v. (60) 



Combining this equation for u with the result for Sv from the appendix (184) we get the 
following, still exact, matrix equation for the perturbations: 



(So £> 2 ) = o. 



(61) 



Here £>o, ©i, and 5^2 contain only terms up to order zero, one, and two in the slow-roll 
functions, respectively. With our choice of the slow-roll functions no higher orders occur. 
They are given by: 

00 -\ o T% + mV T ,-A-B.{4/,4> , ) + a 2 (M?) J - L J' 1 } 

_ ( -H 2 {2~e + f,\\) \ 

^'"^e-VVA (r/lll + e 1 (f) ± )t)(2HP r/ + 3W 2 )y'' [M) 

<°*= n [ f^+^l + W + eift • (64) 



From these equations one can draw the important conclusion that the redefined gravi- 
tational potential u decouples from the perpendicular components of the field <5v up to and 
including first order in slow roll. The resulting equation for u may be written in Fourier 
components (^) as 



3" 



«£+^-tJ*-°- with " S ^T = TO (65) 

From the second expression for 6 and ([54]) it follows that 

J = -H (l + e + r?H) , j = (2e + t?" + 2(7?" ) 2 - (^) 2 - f 11 ) . (66) 
We see that we got back equation (p!9|), which was derived in 1£]. However, here we have 



proved that this equation is not only valid in the single field case, but also in the multiple field 
case, up to and including first order in slow roll, i.e. neglecting the term that is proportional 
to Vli) 1 - . This is precisely one of the combinations of slow-roll functions that is small if 
slow roll is valid. Notice that even if slow roll is no longer valid, the specific combination 
\/e f] 1 - may still be small. This means that slow-roll inflation ends because one of the other 
combinations of slow-roll functions becomes large. We treat an example of this later on in 



section 5.1.2: a quadratic potential where the mass difference between the lightest mass and 
the other masses is large enough. On the other hand, we see that the equation for the field 
perturbations 5v already depends on the solution for u at lower order, although compared 
with other terms this coupling term will become smaller during inflation because it does 
not contain TC 2 . We defined u and 8v in terms of $ and Scf>' L in equation ( |57[ ) with the 
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same powers of e, so that the decoupling holds to the same order for and Scf)^. In the 
literature, e.g. [^4], |2^, ||, this leading order perturbation in u is called the adiabatic mode, 
while the perturbations associated with the perpendicular field components are related to 
the so-called isocurvature modes. In this paper we restrict ourselves to the adiabatic mode, 
but in a subsequent paper we plan to investigate the perpendicular field equations and the 
isocurvature modes. 

The decoupled equation for u cannot be solved analytically in general. However, we can 
solve equation (^35|) for u exactly in the two limits of large and small k, as was observed in 
[fl9f1 . The two linearly independent solutions are denoted by ek and its complex conjugate 
e£. In the small wavelength limit the solution is 



ik{rt-rt k ) 



for 



(67) 



where the time % is defined as the time when k is equal to \Q'L/Qk\- (We use the non-bold 
subscript k to indicate that a quantity is evaluated at r/ = The normalization of the 
two independent solutions is such that their Wronskian satisfies 



In the limit of large wavelengths the solution of (|65|) is 



e k (r?) = C k 6( V ) + D k 6( V ) 



dr/ 



for k < 



which can be rewritten in terms of comoving time as 



e k (t) 



H 
4£l 



+ — ^- D k ( 1 - — f air) di\ . 



(68) 



(69) 



(70) 



Here we have used the definition of 6, the relation e = (1/H)' , and integration by parts. 

Simply joining the two limits in a continuously differentiable way at rj = rjk determines 
the integration constants: 



and 



IX 



tk-f 



(71) 



This joint solution gives a good approximation of the true solution of (pq) if slow roll is 
valid, because then the transition region is small, as we now show. First, we define the 
transition region as that region where the terms k 2 and 6"/ 6 in the equation of motion 
for u are of the same order, i.e. within a factor a 2 of each other with a 2 ~ 10: 



a 



-2i.2 



k Z < 

\0"/9\ < a 2 k 2 . We want to know what the size of this transition region is in terms of r/, 
compared to the characteristic time scale 27r/fe of the solution (]67|). To lowest order in slow 
roll \9"/9\ = a 2 Hl{2e k + ??J) and k 2 = \0%/0 k \ = a 2 k H 2 (2e k + f?J). Hence if we define r?_ as 
the beginning and ij + as the end time of the transition region, we have 



a 



±1 



(Q"/Q) 

(9"/9) 



1/2 



'Ik 



a(v±) 



(72) 
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By integrating the approximation H.(rj) = a'{rj)/a{rj) ~ a{r])Hk around the transition time 
we find that 1 — ak/a(rj) = akH^iji — %). This leads to the following leading order 
expression for the ratio of the duration of the transition and the period of the oscillation: 



77+ -ri- 



le a 



a 



2-n/k a-kHk 2ir 



2efc + % 



a — a 



2tt 



(73) 



Hence if the slow-roll functions are all individually smaller than 0.01, which is a reasonable 
assumption as we will show in the examples in section || then the transition region does 
not last longer than approximately a tenth of an oscillation period. 

We can also discuss the accuracy of our treatment of the transition between the small 
and large wavelength limits in a different way. As we discuss in the next section, the 
important quantity is |-Dk| 2; 



IAcI 




2 h 



(74) 



where in the last step we have included only the leading order term in slow roll. We estimate 



the maximum error in \Dk\ 
matching at and at 77+: 

I Ad 2 - \D k \l 



as the difference between the values of \D^\ determined by 



(\D^)UV--V+) 



\D k \l 



2 [2e k + rj« 



[a 



(75) 



to leading order in slow roll, where we made use of (54) and (JT3[) . Prom this we conclude 
that we can indeed only give the leading order term for |-Dk| 2 i n (|74|), since corrections at 
the next order are expected for a more accurate treatment of the transition region. For the 
single field case an expression for |-Dk| 2 that is accurate up to and including next-to-leading 
order terms was obtained in [27, 17]. 

We finish this section with an argument why the difference between 77/% and rjn does not 
matter in the expression for |-Dk| 2 to leading order. Here tjh is the time of horizon crossing, 
defined by k 2 = TL 2 , which is conventionally used in the literature. On the other hand, rjk is 
the time when the solution of the differential equation (|65| ) changes its behaviour, defined 
by k 2 = \0"/6\ = 7i. 2 (2e + fj^ + ...), which we use to compute |-Dk| 2 - The difference between 
these two expressions for k 2 can be quite large during slow roll. However, just as in the 
previous paragraph, the relevant quantity is 



lAc 12 



H 



ml 



Vk) 



IAcI 2 , 



2V 2e ~fc + € 



(76) 



to lowest order in slow roll. Here we used that an = a^\j2ek + as follows from the 
definitions above, and inserted this into the expression in the text above (^) to calculate 
Vh ~ Vk- We see that corrections to |-Dk| 2 because of the difference between rjk and r\n are 
of higher order in slow roll, if we take the slow-roll functions to be small individually. 
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4 Quantum correlation function of the Newtonian potential 



The quantum correlation function (<3?(x, i])<&(x + r, 77)} during inflation with multiple scalar 
fields is the central object of study in this section. In particular, we obtain an expression 
for the correlation function at the end of inflation. Before the actual calculation can be 
performed various questions have to be addressed. First of all, can one simply quantize 
the Newtonian potential? Another question is which quantum state is appropriate for the 
computation of the correlator. 

In equation (|37|) of section |3.4| we have introduced the variable u and shown that it 
decouples from the perpendicular scalar field fluctations 5v if corrections of the order of 
^rj 1 - are neglected, and that its equation of motion is then given by (|65|). However, the 
Newtonian potential <I> does not correspond to a physical degree of freedom within the 
metric; only the graviton states represent physical degrees of freedom and can be quantized 
in an on-shell quantization procedure. (Alternatively, one can use BRS quantization to avoid 
making explicit gauge choices [ffl, B|.) The scalar field perturbations 5(f) on the other hand 
are physical degrees of freedom, hence they should be quantized. But we are not interested 
in all these scalar perturbations: only those that are directly related to the Newtonian 
potential to this order in slow roll. The relevant multiple field generalization of the variable 



v introduced in [19] is 



0' 

—u 



(77) 



where we have used (^IJ) together with the definitions (^?j) and (|65|). A slightly different 
form of this variable (without the factor a/ n) is sometimes refered to as Mukhanov- Sasaki 
variable [18, p5j| . The equation of motion for v can be found using the exact equation of 
motion for u, i.e. (pO) combined with (|66|), and the expressions for the derivatives of H 



and the slow-roll functions (| 



Av 



(vgr , 
1/6 • 



2V2H 2 



H((3 + e + i] ll )fi + n -Sv + fiVnSv 



(78) 



The definition of v also includes a term with $, which ensures that the equation of motion 
for v can be written in terms of v and Sv only. (If we were not in the situation where \& = <£, 
the definition of v would contain ^ instead of <!>. This definition is automatically gauge 
invariant, as can be seen from (|2|), and therefore it is guaranteed that no non-physical 
degrees of freedom are quantized.) The scale factor a is introduced to remove the first 
derivative term in the equation of motion. The result is that at the beginning of inflation, 
when k 2 S> \(1 / 9)" / (1 / 8)\ as we shall discuss below, the left-hand side of the equation for 
v is simply the equation of the harmonic oscillator. On the right-hand side we see that v 
decouples from 5v up to the same order in slow roll as u does, provided that £ is small 
as well. Hence we know how to quantize v. Moreover, in this limit of large k the equation 
of motion for v is equal to the one for u (pq), up to terms of order ^/Ifj in slow roll. 
This means that the quantum operator u can be expanded in terms of the same creation 
and annihilation operators as v, so that once we have determined the normalization by 
quantizing v, we can simply return to the equation for u to determine the time evolution 
and compute the correlator at later times. 
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In the approximation where v decouples from <5v, the action for v given by 



S = I d 3 xdT/ k 2 



i . ,. 2 i / (l/ey 

(v') 2 + -v(A + 



2 K ' 2 V 1/0 



(79) 



gives rise to the equation of motion ( |78|) . We continue by canonically quantizing the variable 
v at the beginning of inflation rj = rfc. In the quantum theory v and its canonical momentum 
■ft = k?v' satisfy the commutation relation 

[v(x, r]),-fr(y, rj)] = iS(x - y). (80) 

At the beginning of inflation the field v can be expanded as follows: 

* (x ' v) = I (^7E K (??)e " x hl + ^ (r?)£ikx Sk l ' (81) 

where the time independent creation and annihilation operators 6^ and 6 k obey the com- 
mutation relation 

[b k ,b\]=S(k-l). (82) 

In order for this expansion to be valid the dimensionless mode functions v k (r]) and f£(r/) 
have to satisfy the Wronskian condition 

W(vt(r,), vkfo)) = <(77)<(r?) - ttffaKfa) = i 2fc, (83) 

where the v k only depend on the length k of k. It can be checked that the equation of 
motion implies that (v^, v k ) = 0> so that the quantization procedure is consistent with 
the equation of motion; this is guaranteed by using canonical quantization. At the initial 
stages of inflation the scales k 2 that are observable in the CMBR today are much, much 
larger than the horizon: k 2 7i 2 . We assume that this implies that k 2 is also very large 
compared to the other relevant quantities \9"/6\, (9'/9) 2 and \(l/9)" /(l/9)\, since they are 
all proportional to H. 2 , up to factors of order unity or slow roll. Therefore, when inflation 
starts the variable v can be expanded in terms of the independent solutions e k {rj) and e^(r/) 
defined in ([37]). As explained below ([F|), we can expand u in the same way. Introducing 
the notation 

6 = (J , B k = f|j , E k = (e k 4) E k = ikE k @, (84) 

we may write v and u as 

/d 3 k f H 3 k 

<S)i7vlK £kflke " fc9Bk ' !i(x '•' , = / T^ EkU ^' M ^ (85) 



Combining the Wronskian condition for e k (|68|), written as VF(.Ek) = E k QE^ = i2k, with 
the Wronskian condition for v k (||), W(E k R k ) = E^R^QRIeI = i2k, we see that R k is 
an element of the non-compact unitary group U(l,l) defined by 

R k QR[ = 9. (86) 
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In addition the Hermiticity of v implies that is an element of SU(1, 1). The matrix i? k 
represents a Bogolubov transformation HQ. Using the relation ( |77| ) between v and u, ( |84| ) 
for E' k , and the fact that k 3> \0' /9\, the expansion in quantum mechanical operators leads 
to 

u *=mk eR *- (87) 

We now derive a general compact expression for the expectation value of the correlator 
(u(x, rf)vJ(pc+r, rj)) p computed in an arbitary state that is represented by the density matrix 
p. We assume that expectation values (biJ)\) p = 0, etc., if k 7^ 1. Then 

(bl) P i + WpJ ( } ' 

with the number operator iV k = b^b^. The two-point correlator becomes 

(u(x)uUx + r))- - I Eu Uu e~ ikxCT I ^ p ^) I e ik(x+r),T U ] E ] 

{ [ )U { + ))p ~ J (2tt)3 ^ Uk [ (hl) p 1 + (N k ) p ) 6 U * ^ 

We illustrate this expression with two examples. The simplest state to consider is a 
conformal vacuum state density matrix po = |0)(0| where for all k the state |0) is annihilated 
by 6 k . The correlator then takes the form 

/H 3 k 
{A7rk) 3 K 2 1(^)22 e£ - (i? k )i2 e k | 2 e- kr . (90) 

This shows that i? k represents a Bogolubov transformation. In other words, i? k measures 
the alignment of the expansion of the field u in terms of creation and annihilation operators 
(&£ and S k ) with respect to the Lorentz conformal vacuum. In the Lorentz conformal vacuum 
all matrices i? k are equivalent to the identity, so that the field expansion takes the familiar 
Lorentz-invariant form 

" (X ' V) = J (I^k { ei(fc " kX) ^ + e ~ i(fo " kX) ' k l ■ (91) 

Here we used the definition of the functions e k given in ( |67|) for modes with k 2 3> \0"/6\. 
The canonical Hamiltonian H associated with the action (|79j) is given by 



k r 



H = ! d 'k- 

where we wrote 



sinh(2tf) ( e i2S b{bl k + h.c. ) + cosh(2tf) (2iV k + 1) , (92) 



_ /coshtfe^+ 5 ) sinh^e^- 15 ) \ , , 

k ~ Vsinh # e"^-*) cosh e^^ 5 ) J ' ^ 

which is a representation of the most general form for an element of SU(1, 1). Here 1?, if 
and (5 are real, and in general depend on k. We see that for $ = there is no pair creation 
of particles and anti-particles (the b\P_^ term drops out), and cosh(2i9) takes its minimum 
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value at $ = 0. Hence choosing the Lorentz alignment $ = seems to lead to a configuration 
with the usual vacuum properties. 

Our second example is a thermal state of temperature 1//3 in Planckian units charac- 
terized by 

pp oc e~ p6 , H = J d 3 kA;iV k , (N k ) P0 = ^37- (94) 

Here we chose the Lorentz alignment and neglected the (infinite) zero-point energy, which 
is irrelevant for the definition of the density matrix. In the thermal state we have that 
(®L)pp = 0, hence we find for the thermal correlator 

(«(x,77)u(x + r,r,))^ = J ^^ \e k (r,)\ 2 (l + 2(N k ) p ^ e~ ikr , (95) 

where the occupation number {N k ) Pg is given in (|9"4|). 

Next we argue why taking the vacuum state |0) at the beginning of inflation is a rea- 
sonable assumption for the calculation of the density perturbations that we can observe in 
the CMBR today. Even though perturbations in the CMBR have long wavelengths now, 
they had very short wavelengths before they went through the horizon during inflation. 
Therefore, their scale k at the beginning of inflation at tj is much larger than the Planck 
scale. It seems a reasonable assumption that modes with momenta very much larger than 
the Planck scale are not excited at ti, so that for these modes the vacuum state is a good 
assumption|] 

This argument becomes more convincing if we put in some numbers. Assume that 
inflation starts around the Planck time, ti = tp, when the horizon is naturally of the order 
of the Planck scale: H ~ k^ 1 . Of course we do not know the exact quantum state at 
the Planck time, but as a first indication we take a thermal state with a temperature of 
the order of the Planck energy: (3 ~ 1.0 The conformal Hubble parameter is given by 
H = Ka>ie N H ~ cne N with N = J. Hdt, where in the last step we made the approximation 
of a constant H during (the beginning of) inflation. Now we evaluate Ti at the time when 
the scale k goes through the horizon, k = 7i. Suppose that inflation started just before this 
happens, say N ~ 2, and that ~ 1 (the universe should at least be a Planck length at the 
Planck time, and a larger value only makes the argument stronger). Even then we find that 
the thermal correction is already very small: 2{N k ) pg ~ 10 -3 . In physically more realistic 
situations the number of e-folds iV can easily be of the order of 100, so that this thermal 
effect is completely negligible because exp(explOO)) gives a huge suppression. 

As we are making a claim concerning physics at a time when we do not have any solid 
theoretical or experimental data, we have to make sure that our statement does not depend 
strongly on the precise choice of the initial time ij of inflation. However, because \9"/9\ <C 
k 2 the field equation (j6lj) of u is a plain wave equation, it follows that the Bogolubov 
transformation between the Fock spaces at different initial times can be neglected. This 
implies that in our statement there is no fine-tuning problem associated with the starting 
time of inflation. 

5 There could be a problem with this approach, because our knowledge of physics beyond the Planck scale 
is extremely poor. In particular, the dispersion relation w(k) = k that we used implicitly might not be valid 
for large k: there might be a cut-off for large momenta. For a discussion of this trans-Planckian problem 
and possible cosmological consequences see |la |. 

6 Other non- vacuum initial states and possible observable effects have been investigated in Jig], 
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In the remainder of this section we compute the vacuum correlator for the gravitational 
potential <3?, derive expressions for the amplitude and the slope of the density perturbation 
spectrum, and compare our results with the literature. Using the Lorentz alignment for the 



vacuum, the vacuum correlation function for <I> = k \4>\u can be calculated from (9C) and 
®: 

<*(x, t)*(x + r, *))„ = \j ^3 | A<| 2 (l - f jf a(r) dr) e"**, (96) 

where we have used that for large time t >tj. the first term in (|70|) can be neglected because 
of the 1/a suppression. The only model dependence resides in the norm of the coefficient 
|Ac| 2 , given in (fri|) . 

Using the fact that |-Dk| 2 only depends on the length of k, we can perform the integration 
over the angles and obtain: 

($(x,t)d>(x + r,t)) = ^ O °^^|^)| 2 , (97) 



where 

l^(*)| 2 = 4^lAc| 2 (l-f /V)dr) . (98) 
To obtain an estimate of the size of this quantity, we first observe that 

1- — / aM6r=(- [ o(t) dr) . (99) 
a J tk \aJt h ) 

Because H(t) and a(t) are positive functions, we see immediately that the left-hand side has 
to be smaller than or equal to one for t > tk- On the other hand, as long as a(t) is positive, 
and does not grow faster than exponentially, the term within brackets on the right-hand 
side will be a non-decreasing function of t, so that its time derivative is non-negative. Hence 
for all cases of interest 

0< (- f adr\ < 1 \5Ut)\ 2 < ^il^kl 2 - (100) 



a 



i, 



As mentioned in the introduction, we extrapolate our result right to the time of recom- 
bination for the purpose of comparing it with the results in the literature. In 2£] it is shown 



that this is justified for adiabatic perturbations, which are the only ones we consider here. 
However, a full discussion of the validity of this extrapolation is beyond the scope of this 
paper. At the time of recombination during matter domination a(t) oc i 2 / 3 , which leads to 

l<5 k | 2 = — AtIAcI 2 - (101) 

To leading order in slow roll this is exactly the same result as that obtained in [10], if one 
takes into account that the 8"^ defined in that paper equals ||5k| 2 - The calculation of the 



slope of the spectrum is analogous to the one presented in |1C| and gives to leading order 
in slow roll: 

dln|J k | 2 I, . N 

n-l= = -4e fc - 24 = -6e fc + 2 Vk , (102) 
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where in the last step we took the single field limit and switched to the conventional slow- 
roll parameters fl5q) for the purpose of comparison. We see that in that case the result is 
identical to the one in Comparing with the result in p^] we see, after rewriting it in 
terms of our slow-roll functions, that all corrections are indeed of higher order. 
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5 Slow roll with multiple scalar fields 



5.1 Slow roll on a flat manifold 

In this section we consider some examples of slow-roll inflation with scalar fields living on a 
flat manifold. Since all flat manifolds are locally isomorphic to a subset of W N , we assume 
that the N scalar fields live in the M. N themselves. In particular, we use the standard basis 
for M> . The (zeroth order) slow-roll equation of motion and Friedmann equation for the 
background quantities are given by 

= --jLa T v / Vm H = ^y/Vtf). (103) 



We make use of the hat to indicate a unit vector: = 0/0, with (f> = y cf> T cf> the length 
of the vector 0. In the first subsection we consider a quadratic potential where all scalar 
field components have equal masses, while in the second subsection we focus on the more 
complicated case of a quadratic potential with an arbitrary mass matrix. 



5.1.1 Scalar fields with identical masses on a flat manifold 

In this example all masses are assumed to be equal to re m, so that the mass matrix 
is proportional to the identity matrix and the potential reads V = ^K~ 2 m 2 (j) 2 . The mass 
parameter m is given in units of the reduced Planck mass k . The slow-roll equation of 
motion for the background fields simplifies to 

■ * A / 2 777, a / 2 771 X 

= 00 + 00 = -^- — 0=-W-_ and = 0. (104) 

Here we have used the fact that an d 4> are perpendicular, as can be seen by differentiating 
the relation = 1. This means that the direction of is fixed in time; only its magnitude 
changes. The scalar equation can of course be solved easily, and we obtain 

m = (' " h) *<>■ with ( »=vf^' (105) 

where we used the initial condition 0(0) = O . Here too is the time when = if slow roll 
would be valid until the end of inflation. 

Using this solution, we calculate the Hubble parameter H and the number of e-folds N: 

i / \ t / \ 2 

m ± m(j)o ( 1 t \ A7/jA l mjj t\t 1 2l 2 / 1 t 



E ®-7i*-ifV-iz)- m-l^-^-^i--) , doe) 



where = jK 2 ^. Next we calculate the slow-roll functions. Since 0(t) is linear in time, 
7/H and t) 1 - are zero to this order in slow roll. This implies that to this order the decoupling 
of the Newtonian potential from the scalar field perturbations is exact, see (pi]). For e we 
find 
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Clearly, e becomes infinite when t — > t^, which is in contradiction with the bound e < 3 



derived in subsection 13^. But of course slow roll has certainly stopped when t > t% = 
too — V^k/itl, because then e > 1, so that results obtained from equations valid only within 
slow roll cannot be trusted. Notice that we do not have a slow-roll period at all if <po < V2/k, 
since then t\ < 0, so that e is never smaller than 1. We regard the quantity as the 
number of e-folds at the end of inflation. This might seem questionable, as t\ is a better 
candidate for the end time of (slow-roll) inflation than too. However, the difference in the 
number of e-folds using too and t\ is small: — N(t\) = i, so that we can safely use iVoo 
as a good approximation for the number of e-folds at the end of inflation. 

We finish by calculating the model dependent factor | | 2 , which we need to compute 
the correlation function of $ (see section . To this end we must determine the time tk when 



a mode function changes its behaviour (see section ^4). We are especially interested 
in those scales that are observable in the CMBR. As mentioned in the introducion of this 
paper, the corresponding mode functions change behaviour in a small interval about 60 
e-folds before the end of inflation. Hence we consider iV& ~ 60 to be a fixed quantity, and 
determine tk by means of the definition = — N(tk). We find 



£ =i -v£ * h =^k' (108) 

so that slow roll is still a good approximation at time tk- To leading order in slow roll we 
obtain the following expression for |-Dk| 2 : 

W = ^= 2 -m 2 Nl (109) 

5.1.2 Scalar fields with a quadratic potential on a flat manifold 

Now we consider a more general symmetric mass matrix k ^xx\ in the potential. It does not 
necessarily have to be diagonalized, but because it is symmetric we can always bring it in 
diagonal form. As a further assumption we take all eigenvalues to be positive, otherwise 
the potential would not be bounded from below. The potential is denoted by V2 and given 
by 

V 2 = V 2 T m 2 0, (110) 



so that the slow-roll equation of motion ( |103[) reduces to 

2 1 

3 k 2 y/(f) T m 2 (f) 



The solution of this vector equation can be written in terms of one dimensionless scalar 
function ijj(t) as 



0(t) = fH m ^«(£ . (112) 

Here </> = 0(0) is the initial starting point of the field 0, which implies that ip(0) = 0. In 
other words, we have determined the trajectory that the field <\> follows through field space 
starting from point O - 
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An important role in our further analyses is played by the functions F n , defined by 

Fn = ^ -2 — , (H3) 

n 

with 4>q the length of 4> : 4>q = </>o <^o- The functions F n (ip) are positive and monotonously 
decreasing for all ip, tending to zero in the limit ij) — » oo, because we have assumed that 
all mass eigenvalues are positive. Using these definitions we see that the function ip(t) is 
determined by the following equations: 

V 3/^00 y/Fi(ll>) 

Notice that tp is always non- negative. The functions F n do not depend on the length of (f> , 
only on its direction, as can be seen from the definition ( |113 ). This implies that the only 
dependence on the length <po in ( p.l4j ) is in the factor 2/(k 2 0o), so that it can be absorbed 
by a redefinition of the time variable only. 

Next we discuss some additional properties of the functions F n . The definition of F n 
can also be written as 

F n = 72 > (ll 5 ) 

for any integer — n < p < n. Using the Green-Schwarz inequality (A T B) 2 < (A T A)(B r B) 
for arbitrary vectors A and B, we obtain 

F^ < F n _ p F n+P . (116) 
From the definition of the F n we also see that 

^nW) = -fn+iW- (117) 
We can express many important quantities in the functions F n . The Friedmann equation 



( 103 ) for H simplifies to 



H = ^= y/fW)- (H8) 



Formally integrating equation (114) for ip, we find how long it takes to go from ip = to ifi 
in slow roll: 

1> 



%) = a/o k Vo/ d^'v^j. (119) 



The number of e-folds N = J Hdt can be interpreted as a function of ip given by 



Nty) = ^ I atfFiW) = Noo (1 - FqM) , (120) 



where we combined (11?) and (|119| ), and used (|117| ) to perform the integration. The number 



of e-folds in the limit ip — > oo is given by Noo = 4k 2 0q, which can be interpreted as an 
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upper limit for the total number of e-folds during slow-roll inflation. Finally, the slow-roll 
functions (|55h can now be written as 



2 F 2 „__ 2 F3F1 - F$ ._ 2 VF 4 F 2 - jg 

'/ — ..9.J.2 rfrrn ' " ..5 12 77-77. • 



K 2 </> 2 F 2 ' ' K 2 2 F 2 F2 > / K 2 2 ^ 



Using the Green-Schwarz inequality (116) we see that fj" is always negative, while fj is 
real, as it should be. Observe that if m is proportional to the identity, the inequality is 
saturated and 77" and fj 1 - are zero. This is in agreement with the results of the previous 
subsection. Since the functions F n (tp) are independent of <pQ, this dependence enters only 
in the prefactors of the expressions for t(ifj), N(ip) and the slow-roll functions. 

Before going on to discuss estimates for the functions F n , we need to introduce some 
additional notation. We define a semi-positive definite matrix norm: 

l|A|| 2 = ^, (122) 

for any arbitrary N x iV-matrix A. The reason that || • || does not define a regular norm 
is that ||A|| 2 = does not imply that A = 0; we can only infer that A0 O = 0. Indeed, if 
A has determinant zero and cf> is one of A's zero modes, A</> = is satisfied without A 
being the zero matrix. With this norm the definition of F n (ip) can also be written as 

F n {$) = ||m n e-2 m ^|| 2 . (123) 

We order the eigenvalues of m 2 from smallest to largest, m\ < m| < . . . < m|. Here we look 
only at distinct eigenvalues, so that £ is smaller than N if there are degenerate eigenvalues. 
The projection operator E n projects on the eigenspace with eigenvalue m 2 . These operators 
are mutually orthogonal and sum to the identity: E ra = 1. The norm of these projection 
operators satisfies ||E n || < 1. 

Above we have been able to write all kinds of important quantities for the slow-roll 
period in terms of the functions F n (ip). But these functions are rather complicated as they 
depend both on an (exponentiated) mass matrix m 2 and on the direction of the initial 
vector cf) . Since in the CMBR we cannot see further back than about the last 60 e-folds 
of inflation, while the total amount of inflation is generally much larger, we are often only 
interested in the asymptotic behaviour of the functions F n {ip) for large ip. Below we will 
show that the asymptotic behaviour is indeed a good approximation for the time interval 
during inflation that is (indirectly) observable through the CMBR, but first we concentrate 
on the asymptotic expressions themselves. 

As can be seen from the definition of F n in ( |113| ), in the limit ip — > 00 the smallest mass 
eigenvalue will start to dominate. We denote the smallest eigenvalue hy p, p = mi, while 
the ratio of the next-to-smallest and smallest masses squared is called p: p = m 2 /m 2 > 1. 
Furthermore, the operator E = Ei projects on the eigenspace of the smallest eigenvalue^] 

Here we assumed that ||Ei|| 7^ 0. If ||Ei|| is zero, it means that </> has no component in the sub spac e 



7 



corresponding with this eigenvalue. As can be seen from the differential equation for this quadratic case (111), 
this means that <j> will never obtain a component in the directions corresponding to this subspace. Hence 
we remove these directions from the problem, consider ^ to be a vector of appropriate (lower) dimension, 
and take mi to be the smallest remaining eigenvalue, etc. 
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and we define x = l|E2|| 2 /l|Ei|| 2 - Using these definitions, we find the following asymptotic 
behaviour for the functions F n {ip) in the limit ip — > oo: 

F n -> \\B\\ 2 i2 2n e-^ (l + xp n e- {p - 1] ^) -> | |E| | V^e"^, (124) 

where the first limit contains both leading and next-to-leading order terms, while the second 
contains only the leading order term. Both these asymptotic expressions for F n are needed 
to obtain the non-vanishing leading order behaviour of ratios and differences of ratios of the 
functions F n : 

f - » 2ip - q \ _ III- _> M 2 p n-l (p _ I) 2 xe -(P-1)MV. (125) 

Using these expressions we find the asymptotic behaviour for the Hubble parameter ( [IT 
and the number of e- folds ( |120| ), 



ff^^^He-W, ^^^(l-llElpe-^). (126) 



The asymptotic behaviour of the slow-roll functions (|121| ) is given by 



K 2 0§||E|| 2 ' ' k 2 <$||E|| 2 



2 V* ^(p-^e-l^- 3 )^. (127) 



k 2 2 ||E|| 2 

Notice that 77" goes to zero for p > 2, while for p < 2 it diverges. The same holds true for 
f] ± , but there the critical value is p = 3. Since p > 1 by definition, the slow-roll function 
e always grows faster than 77" and rj 1 - in the limit ^ — > 00. As we discussed in section 3.3 , 
however, it is really the combinations y/tfjW and y/if/ 1 - that determine whether slow roll is 
valid or not. Their critial values are p = 5/2 and p = 4, respectively. 

We finish this subsection by calculating the model dependent factor |-Dk| 2 . Following 
the same steps as in the previous subsection we find that = tp{tk) and = e(ifc) are 
given by 

'-^-Pii^ *-w t - < 128 » 

So as long as S> our assumption of using the asymptotic behaviour for ip —* cxd 
at time is very good. Moreover, since in this limit e is the largest of the three slow- 
roll functions, we see that the slow-roll approximation is also valid. Hence we can use the 
leading order slow-roll estimate of fl74]) for |-Dk| 2; 

K 2 H 2 9 

\D^ = ^ = {p 2 Nl (129) 
2 ek o 



where, apart from the previous two expressions, we also used (|126| ) for H. This result agrees 
with (|109) for identical masses. 
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5.2 Slow roll on a curved manifold 



Now we turn to the slow-roll behaviour of scalar fields that parameterize a curved manifold 
that is isotropic around a point. We start with setting up the general framework, which we 
clarify by examples and expand upon in the special cases discussed in the next subsections. 
Consider an TV-dimensional manifold with coordinates and metric G(0) given by 

G(0) = 5 (0) (l N + J^j^ Q) , (130) 
with g(4>) 7^ and A(0) ^ 1. The matrix Q is the projection operator defined by 

Q = e eQ , with e = ^. (131) 



Here <fi = \J 4> T 4> repr esents t he coordinate length of the vector 0, which should not be 
confused with |0| = y T G0. By taking this form for the metric we have of course re- 
stricted ourselves to manifolds that are isotropic around a point, but it covers some general, 
interesting cases, e.g. the sphere with either embedding (see section [5.2,2 ) or stereographical 



coordinates. Also, the equations of motion in a central potential obtained with this metric 
are identical to those obtained in the case of a more general metric, as we will explain in 



section 5.2.1. The inverse of this metric and the determinant are given by 



1 a N 
G 1 = - (1 N - AQ) , detG = -^— . (132) 

g 1 - A 

For the determinant we used the relation lndet G = trlnG and the fact that tr Q = 1. 

Inserting our special choice for the metric into the (pure) slow-roll equation of motion 
for gives: 

= —?L G~ l d T ^VW) = ~— (d T W - ^(dW 0)0) . (133) 
V3k v3Kff V <P J 



Notice that dyV = 4> a d a \/V is a scalar. In general this vector equation can be hard 
to solve, but in practice we often have some information from the corresponding flat case 
that we can use. In particular, we can often determine the trajectories that the scalar fields 
follow through the flat field space. On the other hand, it is much harder to calculate exactly 
how the scalar fields move along these trajectories as a function of time, but this still means 
that we have reduced the system of N differential equations for to a single one that gives 
the velocity along the trajectories. In other words, the trajectories of the slow-roll equation 
of motion for the flat case can be written as 

flat (t) = T(^(i),0 o ), with T(O,0 o ) = 0o and V(*o) = (134) 

(with T a known function), where the function t/j(t) has to satisfy the differential equation 

V* = - I (flat case). (135) 



An example of this was given in ( |112| ) and (114) for the case of a quadratic potential. 
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This flat solution can be generalized to curved manifolds with a metric of the form 
introduced above by defining 

0curved(*) = TW), 4> ). (136) 

Here T is the same function as above, while the differential equation ( |135| ) for if) is slightly 
modified to 

*-»i-^t (curved case). ,137, 

By inserting our ansatz for the solution into the equation of motion we find that the factor 
s(ip) has to satisfy 

" M = A^1^^5^ and ,(0) = 1. (138) 
a T T T dy/VTrf 

We give examples of this general method in the following subsections. 

Next we discuss the definition and evaluation of the slow-roll functions e, fj" and j)^. To 
this end we define functions C n (V) as follows: 

Ci(V)=k C n (V) = k — — ^~~PT~ > n>2. (139) 

n n 

The functions C n (V) are more than simply the curved generalization of the functions F n : the 
C n are defined for an arbitrary potential, while in the definition of the F n we have assumed 
a quadratic potential and made use of the fact that we can determine the trajectories of 
the fields in that case. Using the Green-Schwarz inequality we can derive the following 
inequalities for positive integers n,p with < p < n: 

C 2 < C2pC 2 ( n - p ), (140) 

which follows by writing <f%C n = K 2n Vy(G- 1 V T Vy) ! '- 1 G- 1 (V I 'VFG- 1 )' ! ^ 1 V T F. 
The slow-roll functions (M) can be written as 



which are the same expressions as in the case of a quadratic potential on a flat manifold 
(121), but with the F n replaced by C n (V). The only inequality for the functions C n that is 
directly applicable is for n = 3,p = 1: Cj < C2C4, which implies that the square {fj L ) 2 is 
positive, as it should be. 

On a curved manifold the second order covariant derivative V T Vy contains connection 
terms. We now compute what these terms are in the case of our special metric. It is 
convenient to work out the special combination 

v = (V T VV)G- 1 V T V, (142) 

since that is how the connection enters into the expressions of the slow-roll functions, as 
can be seen by writing 

C 3 = K ^ VK Z v , C 4 = k 8 V ^ V . (143) 

n n 
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This vector v can be split into two vectors v = Vp + vg, given by 

v F = (d T dV)G- 1 d T V, (v G ) a = -F c ab d c V(G- 1 d T V) b = hiG-^&VdcV. (144) 

Here we have used that -T c ab w c {G- 1 w) b = -T cab (G~ l w) b (G~ l w) c = ^(G~ l ) hc a WbW c for 
any vector w. Next we compute the derivative of the inverse metric. Using the definitions 



of the inverse metric G 1 (132) and the projection operator Q (131), we obtain 

Here we have assumed that the metric functions g and A only depend implicitly on cf) via 
the quantity x = 4> 2 /R 2 , with R the characteristic radius of curvature of the manifold. This 
leads to the final result for vg: 



1 



gR? 



9 , xd VG-WV + ^*)(8V*f 



4,--(dV4,)^f. (146) 
9 V 



In the next subsections we work out these expressions for the slow-roll functions in the cases 
of some special potentials. 

5.2.1 Scalar fields on a curved manifold with a central potential 

We consider the case that the potential V c {4>) is a central potential around the origin: it is 
a function of the coordinate length cp only. The first and second order gradients d T V c and 
d T dV c are then given by 

d T V c = V c ^e and d T dV c = V C ^Q + ^V c ^(l - Q) . (147) 

As in our first example of scalar fields with equal masses on a flat manifold in section |5.1.l| , 
we find that the vector slow-roll equation of motion reduces to a scalar equation: 

<t>= — 75 7=~eo => = — 75 7=r- ( 14 8) 

V3k 9 VV C V3k 9 VV C 



Notice that if we take a more general metric 

/ N-l 

K(<t>) 



with a sum of mutually orthogonal projectors Q n , instead of just the one Q = Qo, no 
other terms appear. This is because dV c is pointing in the radial direction eo and eo is an 
eigenvector of the metric G as well as of metric G: Geo = Geo = j^eo. 

Before going on to discuss the slow-roll functions to leading order in slow roll, we note 
that we can also say something about the exact slow-roll function fj^. From the exact 
definition of r) in (^) and the exact equation of motion (|47]) we get 

= avd-e^G-^v = i-A rv^_V _ 

U H^\ 2 9 \H\<f>\J K h K 1 
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where a denotes the angle between the unit vectors eo and ej, or equivalently, between the 
position vector 4> an d the velocity (f>. This angle a is defined in curved field space as 



cos Q 



gi • e 
l e i|l e o| 



1 - A 



e^eo. 



(151) 



Notice that while |ex| = 1, this is not the case for |eo|, as eg was defined as having coordinate 
length equal to one. From this formula we infer that if the field velocity </> is pointing in the 
same direction (up to orientation) as the coordinate vector (j> (i.e. a = 0,ir), the slow-roll 
function rj 1 - vanishes. Notice that if this happens, it holds for all time, as can be seen from 
the exact equation of motion. 

Next we work out the expressions for the slow-roll functions to leading order in slow 
roll, given in (141). As discussed in the previous section, it is useful to first calculate the 
vector v = vp + vq, defined in Ql42j) ff.. We find in the case of a central potential 



1 - A 



Vr._A V, 



g,x <p 
g R 2 



Vc,d> e 



9 R 2 



{V c ,<t>) 2 eo, 



which leads to 

1 - A (V, 



2.9 



ItVr. 



1 - A fV r 



9 



9.x 



k 2 V c g R 2 k 2 V c 



+ 



(152) 



(153) 



Furthermore, fj = 0, as can most easily be understood by realizing that the Green-Schwarz 
inequality is saturated (C| = C2C4) if v and VF are (anti-)parallel. This is in agreement 



with the result (150), since to leading order in slow roll cos 



a 



1. 



To illustrate various aspects of the general discussion above, we now turn to an example: 
a quadratic central potential of scalar fields with identical masses k _1 to, which are the local 
embedding coordinates on an iV-dimensional sphere with radius R. These coordinates </> 
are induced by embedding the sphere in an N + 1 dimensional Euclidean space, so that by 
construction cj) 2 = 4> T cj> < R 2 . At least two of these coordinate systems are needed to cover 
the whole sphere. In the slow roll discussion here, we stay within one coordinate system 
because the quadratic potential is minimal in the origin of this system. The metric, its 
inverse, and the connection are given by 



U + 



R 2 



Q, 



1 



and T a h 



be 



(154) 



Hence g = 1 and A 



' /R < 1 in terms of the general metric (|13 



Inserting the relevant quantities into the slow-roll equation of motion ( |148| ) we find 



2m/ (p 2 
3 ^ { R 2 



(155) 



Solving this equation, multiplying by a constant unit vector, and applying the initial con- 
dition <^>(0) = cj} Q , we get the following answer for the solution of the background equation 
to leading order in slow roll: 



4>(t) = R tanh[ 7 (t 00 - t)\ 0, 



0- 



(156) 
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with 7 = yf and = ^ In ■ The slow-roll functions follow immediately from 



(153): 



R 2 - 4>1 ( sinh(7i c 



W = ^nr^ t = o. (157) 



k 2 R 2 4% V smh [7(^oo -t)]J ' ' K 2 R 2 ' 

Since R is fixed by the model, the slow roll conditions, e and y/tffi small, are satisfied if (fro 
is such that R 2 — (/>q <C R 2 . From the Hubble parameter H = fcyF/3 = mcj)/y/6 we derive 
the expression for the number of e- folds: 

N{t)= / Hdt = N 00 — Incosh^oo-t)], (158) 

Jo 1 

with iV^ = -i« 2 i? 2 ln(l - ^/i? 2 ) > 0, since < < i? 2 . Therefore, the field should 
start close to the equator of the sphere, R 2 — 4>q <C R 2 , to ensure that sufficient inflation 
is obtained. This is compatible with the requirement that the relevant slow roll functions 
(e and y/ifj") are small initially. In the limit R — > oo all results agree with those we found 
in the flat case in section |5.1.1| . Notice that we can also determine the solution for <fi by 
using our knowledge from the flat case and the method described in ( |136j ) ff., but in this 
particular case that is more complicated. 

We finish with the calculation of the expression for the coefficient |-Dk| 2 j following the 
steps outlined in section |5.1.1| . The time and the slow-roll function e at tk are given by 

(2_/V \ 2 1 

* e ~ fc = ^ W 4 M-i ' (159) 

This confirms that the slow-roll approximation is still valid at tk- For the model dependent 
factor |-Dk| 2 we obtain to leading order in slow roll 

|Ac| 2 = |™ 2 f^sinh^V- (160) 



3 V 2 k 2 R 2 , 

In this calculation we have assumed that 77" itself is small, which means that R is larger 
than the Planck radius k . If this is not the case, more terms in ( |74| ) have to be taken 
into account. In the limit that R — > 00 this result agrees with the flat case (|109|). 



5.2.2 Scalar fields with different masses on a curved manifold 

Next we consider the case of the quadratic potential V2, defined in (|11C|) , on a curved 



manifold with metric (130). The slow-roll equation of motion in this situation reads 



= / ( m 2 - 0. (161) 

Since we know the trajectories of <j) m the flat field case (|112| ), we can use the method 
described in ( |136| ) ff. to determine the solution, or at least the trajectories, of <fi in the 
curved field space. The equation for s(ip) can be solved analytically for the sphere with 
embedding coordinates, which was also considered in the previous subsection. 
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Before turning to this example, we first work out the expressions for the functions 
Cl(V%)j-- - , C^V^) in terms of the flat functions F n defined in ( |113| ). These expressions 
can be used to determine the slow-roll functions (|141|). Because we found many properties 



and estimates for the F n in section 5.1.2 , a lot of additional properties are obtained for 
the functions C n (V2) in this way. To find the relation between the C n (V2) and the F n 
it is convenient to define intermediate functions F n that incorporate some of the non-flat 
metric aspects, but not the full covariant derivatives (connection terms). They are defined 
as follows: 



cj) 1 G(G~ 1 m 
ST 



n> 1. 



(162) 



Notice that apart from the metric aspects, the F n also contain an extra factor of s 2 as 
compared to the F n , since = s(ip)e~z m ^</> according to ( |136| ) . The functions F n can 
be expressed in terms of the F n ; for the first four functions F n we find by inserting the 
definition ( |i~32] ) of G 1 



Fx = s'Fr 



F 2 = — (F 2 - AFi) 
9 



F 4 - A ( 2F 3 + ^ ) + 3A Z F 2 - A d iq 



h = -J (*3 - 2AF 2 + A 2 Fx) , 

g 2 



with 



A = A 



F\ 



(163) 



Here we used that ^ T m 2 PQm 2, ?Q . . . m 2r cj> = (Jgg 



as 



F n = r (m 2 G- 1 ) p G(G- 1 m 2 ) ?1 - 



jr )Fos 2 . Writing the functions F n 
, we obtain the same Green-Schwarz inequalities 



as for the C n : 



(Fn) 2 < F2„F<){ 1 



2p^2(n-p) 



(164) 



for integer < p < n. The next (and final) step is to write the functions C n (V2) in terms 
of the F n . It is easy to show that 



C 1 (V 2 ) = F 1 = s 2 F 1 



and 



C 2 (V 2 ) = F 2 = - (F 2 - AFi) . 
9 



(165) 



For the functions C3 and C4 we use the vector v defined in ( |142j ), which in the case of a 
quadratic potential can be written as 



1 



m 2 G W 



A 



H 



-m 



1 0, 



(166) 



where H = g, x %F 2 + AF X ^§ 
and C4 we obtain 



1 



9 9 

. By inserting this into the expressions for C3 



C 4 (V 2 



c 3 (v 2 ) -- 

A- /A\2 
2-F3+ - x 
9 v 5 



9 



1-XH 



Fo 



9 9 
1-XH 



Fi, 



(167) 



9 9 



p 2 p 1 

9 



H 

2fJ 



S 2 Fn 
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Next we discuss the example mentioned above. We consider the case where the manifold 
on which the fields live is a sphere with radius R. We again use embedding coordinates, i.e. 
g = 1 and A = x = (f) 2 /R 2 . Solving equation ( |13§| ) for s we find 



2R 2 



1 



s 2 (V>) 



1 



R 2 



(168) 



The trajectories of <f> and the differential equation for ip ( |137| ) are given by 

^2 2 1 1 



0) 



(169) 



As in the flat case, ?/> is monotonously increasing, starting at zero on t = 0, and reaching oo 
when = 0. The slow-roll functions are again determined from (|141| ), using the expressions 
for the C n derived above. However, since the resulting expressions are quite large, we only 
give e here: 



1 F 2 



e 2 $ s 2 F 2 k 2 R 2 K 2 (f> 2 V R 2 J F 2 k 2 R 2 F{ 



2 F0F2 — F 2 



(170) 



which, using the Green-Schwarz inequality ( |116| ), is seen to be larger than or equal to zero, 
as it should. For the number of e-folds we find 



N 



2 a2 



s 2 F 1 di(j = N 



1 



k 2 R 2 In [ 1 + 



R 2 



(171) 



where is the same as in the previous subsection: = \k 2 R 2 ln(R 2 /(R 2 - </>§)). 

As mentioned above the functions F n satisfy the estimates (124). This means that we 
can give the following estimates for H, e and iV in the limit ip — > 00, using that <j)Q < R: 



H 



<M|E| 



V6 



N^N n 



1 



R 2 



k 2 R 2 In ( 1 + 



|E||V 



1 R 2 - 4>l 2^ 

1 



k 2 ^||E|| 2 R 2 



R 2 



|E||V 



(172) 



Following the same steps as in the previous subsection to determine |-Dk| 2 ; we find for the 
critical value ipk 



1 R 2 -4>1 



IEII 2 



exp 



4A fc 

K 2 i? 2 



1 expl^ 



1 



l E H 2 ex P (^)-l- 



(173) 



Hence, as in the flat case, the assumption of looking at the asymptotic behaviour for tp at 
time tk is a good approximation if iVoo 3> Nk- Finally, the expression for |-Dk| 2 in this limit 
is 



I A. 



k 2 R 2 . 2A fe 
— - — sinn 



K 2 R 2 



(174) 



in agreement with (160) for identical masses. 
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6 Conclusions 



In this paper we have analyzed scalar gravitational perturbations on a Robertson- Walker 
background in the presence of multiple scalar fields that take values on a (geometrically 
non-trivial) field manifold during slow-roll inflation. 

We have modified the definitions of the well-known slow-roll parameters to define slow- 
roll functions in terms of the Hubble parameter, background field velocity and their deriva- 
tives in the case of multiple scalar field inflation. This means that the slow-roll functions 
have become vectors, except for e which is a derivative of the Hubble parameter. Like other 
relevant vectors they are split into a component parallel to the scalar field velocity, and 
components perpendicular to this velocity vector. To define the slow-roll functions we do 
not need to make the assumption that slow roll is valid, but if it is valid one can expand in 
these functions, giving the relative importance of terms in various equations. For example, 
if e, \/! fj 1 - and \/l fj" are small, the background equation of motion for the scalar fields can 
be approximated by the (pure) slow-roll equation. 

We set up the combined system of gravitational and matter perturbations in a way 
analogous to Mukhanov et al. |19U , but including multiple scalar fields and effects of a 
non-trivial field geometry. The component of the scalar field perturbations parallel to the 
background field velocity can be eliminated. The remaining perpendicular components of 
the field perturbations and the gravitational potential are described by coupled differential 
equations (|6l|). However, the gravitational potential decouples from the perpendicular field 
perturbations if effects of the order of v / ^??" L can be neglected; to the same order the 
background equations reduce to the slow-roll equations. 

Since to first order in slow roll the equation for the gravitational potential is equal to 
the single field case, it has the same expressions for the solution as |n| in the small and 
long wavelength limits. Using the slow-roll functions the corrections due to the transition 
region between these two limits can be estimated. It follows that simply joining the two 
solutions together yields a good approximation during the complete inflationary period, in 
the sense that corrections are of higher order in slow roll. The transition between the two 
solutions happens when k 2 = \G" /0\ ~ Ti 2 (2i + f}^). Even though this is unequal to the time 
of horizon crossing k 2 = Ti 2 of a scale k, the corrections to the gravitational potential if one 
would use this time instead are suppressed in slow roll. 

The quantum two-point correlation function of the gravitational potential is related to 
the temperature fluctuations that are observed in the CMBR. The only physical degrees of 
freedom that can be quantized in the system we consider are the scalar field perturbations. 
(Only the graviton states are physical degrees of freedom on the gravitational side. However, 
they are decoupled because the Einstein equations have been linearized.) Although the 
quantization at the beginning of inflation involves the scalar field perturbations, after the 
Fock space has been constructed the time evolution of the correlator can be calculated as 
if it is a classical quantity. Choosing the vacuum as the initial state of inflation is a good 
assumption for the scales that are observable in the CMBR, since their momenta in the 
initial stages of inflation are much larger than the Planck energy. Taking a thermal state 
with the Planck temperature at the beginning of inflation as a first attempt at improving on 
the assumption of a pure vacuum leads to corrections to the vacuum correlator of the order 
of only 10 -3 or less. Here we assumed that at least a few e-folds of inflation occured between 
the Planck time and the moment when the observable scales went through the horizon. The 
gravitational correlation function contains one directly model dependent factor that can be 
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expressed in terms of the slow-roll functions. 

Finally, we discussed some multiple field examples to illustrate some dynamical aspects 
of the background and compute this model dependent factor. On flat manifolds with central 
or quadratic potentials the trajectories of the fields in field space can be found in terms 
of one function ip(t), so that only a single differential equation remains. This equation 
of motion, as well as many other relevant quantities like the slow-roll functions and the 
Hubble parameter, can be expressed in functions F n of this ip. These functions have many 
useful properties which make it possible to analyze slow-roll phenomena without having to 
explicitly solve the (complicated) equation for -ip. They satisfy inequalities that follow from 
the Green-Schwarz inequality. Their asymptotic behaviour for large ip can be used to analyze 
various important quantities when observable scales go through the horizon, provided that 
there have been many e-folds of inflation before that. Using that the derivative of F n equals 
— F n+ \ makes it possible to integrate the Hubble parameter and obtain the number of e-folds 
as a function of if). 

The generalizations C n of the functions F n to curved field space and to an arbitrary 
potential lead to the same expressions for the slow-roll functions as in the flat case with 
a quadratic potential. For such a potential on a manifold isotropic around the origin we 
expressed the C n in terms of their flat relatives. If the trajectories of the fields in the flat 
field space are known for a given potential, the corresponding trajectories on an isotropic 
manifold are the trajectories in flat space multiplied by a scalar function s of the parameter 
ijj. The sphere with embedding coordinates is a special example of an isotropic manifold. 
The background equations can be solved explicitly as a function of time for a quadratic 
potential with all masses equal. If not all masses are equal, it is still possible to find 
an integrated expression for the number of e-folds in terms of Fq(iI>). For the flat space 
examples we considered, the only possibility to obtain a large total number of e-folds is by 
taking large initial field values. The radius of curvature of a curved manifold is an additional 
parameter that influences the total number of e-folds: for the examples of the sphere this 
number becomes large if the initial field values are of the same order as the radius of the 
sphere. 
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A Derivation of the equation of motion for the perpendicular 
field perturbations 



This appendix is devoted to the derivation of the equation of motion for the perpendicular 
field components Sv = -^rrr Scb . We start by taking the perpendicular projection of 
equation (^) . Using the antisymmetric properties of the Riemann tensor and the fact that 
we contract with two identical vectors cf)' we conclude that P^R(<^', cf>')8cf> = R(0', cf>')8cf)^. 
The perpendicularly projected equation of motion for Scf) can be written as 

P 1 {pi + 2HV n + a 2 M 2 ) + 2a 2 $P ± G~ 1 V T V + P 1 (V 2 + 2HV V ) 8^ 



+ a 2 (M 



2\±± 



A - R(0', 4>') ) 84 >± = 0. (175) 



Here we split 8(f) = 8(f)" + Scf) 1 ' into components parallel and perpendicular to cf)' . Further- 
more, we used that the projection operator commutes with A and employed the notation 
defined in (|36|). The last term can be trivially rewritten, so we now derive expressions for 



the first three terms on the left-hand side of (175) in terms of u and only. 

We consider the first term in (175). This term does not vanish in general, since the 
derivatives and mass matrix generally change the direction of Scf)^ . The covariant derivative 
T> 2 applied to the non-contracted vector cf) in 8cf>^ = ^rm^ cp' is canceled by the mass term. 
This can be seen from the covariantly differentiated background equation of motion (^2|) , 

V 2 cf)' + 2(H'- 2H 2 ) cf)' + a 2 M V = 0, 



(176) 



where we used V rj (G~ V V) = M cf) , if we apply the perpendicular projection operator: 
P^D 2 */)' + P- L (a 2 Nl 2 cf>') = 0. The remaining contribution of the first term of ( |175[) can be 
written as 



P 1 (V 2 + 2HV r , + a 2 M 2 ) Scf^ 



K 2 \cf)' 



K 2 



cf)' -Scf) 



\<t> 



'12 



(177) 



The factor involving cf) ■ Scf) can be rewritten in terms of u using the (Oi)-component of the 
Einstein equation (p|): 



<f>' ■ Scf) 
10 



'12 



2k ( . \cf)'Y 

W\^ U+ W\ U 



nv2i 



u" + 



W 
\<t>'\ 



, (178) 



where we have used that k\4>'\ = Hv2e, which follows from ( p9|) and (|53f) . 

For the second term in fll75|) we use the background field equation (f42|) and obtain 



2a 2 $P- L G~ 1 V T V 



-2$ {V^Y 



K 2 \cf)'\ 



nV2i 



(n' - n 2 )i 



gW) 



>\-L 



(179) 



using another, but equivalent expression for k\4>'\: k\<P'\ = 2(7i 2 — H') / (TlV^e) . Combining 
these expressions and using the equation of motion for u ( |58|) we find that the first two 
terms of (11751) can be written as 



P- 1 (. . .) Scf)W + 2a 2 $P ± G~ 1 V T V = 4 



K 2 \cf)> 



Au 
7iV2e 



(p^y ■ 



(180) 
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Moving to the third term of ( |175| ), we first note that substitution of 8(f)^~ = K jr-^ v 
gives 



.21 Jj\ 



(1-Pll) 



Sv. (181) 



Since P" only gives a non-vanishing contribution if it is applied to T> v dv or D 2 Sv, we 
consider those terms separately and find 



P [V 2 v + 2 l -^-V r ,)Sv 



dv + 2 — ptj; — 'P^v ei. (182) 



10' 



10' 



To obtain this expression, we wrote P" = eicj)'^ /\<fi f \ and used the two relations that are 
obtained by taking the first and second derivative with respect to conformal time of the 
equation </>' ■ Sv = 0: 



' . V^dv = -V n 4>' • 5v, 4> ■ V*8v = -V*4>' • <5v - 2(2^0') • V v 



5v, 



where in the second relation we used V^cf)' = (V^rf)')^ + ^Wl ^' an< ^ ^ e ^ rs ^ re ^ a ^ on - 
Now we can insert the expressions ( |180 ), ( |181|) and ( 182 ) into (|175| ) and find 



? \<P'\' 



+ 4 



1 



nV2i 



Au + 



<t>'\ 



n'-n 2 )+ a 2 (M 2 ) ±x - A - R(0', 0') 



<5v 



/\_L 



I0'l 



+ 



10 I 



(183) 
ei = 0. 



Using (|53| ) we obtain the final form in terms of the slow-roll functions: 



V 1 + 2H[l + ^)V ri + H I [ W + (V ) + £" + a 2 (M 2 ) x± - A - R(0', 0' 



Sv 



+ 



W 2 (3r) + X • <5v + 2?^ • P^v 



ei +4 



W 2 ™ 1 - • <5v + -= An fi^ = 0, (184) 
s/2e i 



which corresponds to the <5v component of the matrix equation (pip. 
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